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Abstract 

It has been observed that the composite fermion (CF) approach tends to overcount the number 
of linearly independent candidate states for fixed sets of quantum numbers [number of particles, 
total angular momentum, and (pseudo)spin if applicable]. That is, CF Slater determinants that 
are orthogonal before projection, may lead to wave functions that are identical, or possess linear 
dependencies, after projection. This has been pointed out both in the context of rotating bosons 
in the lowest Landau level, and for excited bands of the (fermionic) fractional quantum Hall effect. 
We present a systematic approach that enables us to reveal all linear dependencies between bosonic 
compact states in the lowest CF “cyclotron energy” sub-band, and almost all dependencies in higher 
sub-bands, at the level of the CF Slater determinants, i.e. before projection, which implies a major 
computational simplification. Our approach is introduced for so-called simple states of two-species 
rotating bosons, and then generalised to generic compact bosonic states, both one- and two-species. 
Some perspectives also apply to fermionic systems. The identities and linear dependencies we find, 
are analytically exact for “brute force” projection in the disk geometry. 
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I. INTRODUCTION 


Throughout the history of quantum Hall physics, a particularly successful line of research 
has been the construction of explicit trial many-body wave functions, notably the famous 
Laughlin wave function!^^, the phenomenology of composite fermiond^, and more recently var¬ 
ious schemes applying to non-Abelian stated®. The idea is that while not exact [at least not 
in the case of Coulomb interaction], the thus obtained wave functions capture the important 
topological properties of the state at hand. Many of the methods developed in the context of 
the quantum Hall effect have recently been applied to cold atom systems. This is of interest 
thanks to the impressive experimental developments in generating artificial magnetic helds 
in atomic Bose condensates, either by rotatiorl® or other methodd®'. The hope is to be 
able to realize strongly correlated states of the quantum Hall type in cold atom systems. 
This would provide a setting where parameters like disorder and interaction strength are 
controllable and tunable to a much larger degree than for electrons in semiconductor het¬ 
erostructures and might, eventually, be superior to the electronic quantum Hall system for 
e.g. studying topological quantum computing. 

As mentioned above, the composite fermion approach has been a great success in describ¬ 
ing a large number of quantum Hall states, including their fractional excitationd^. Later, 
this formalism was modified to successfully describe rapidly rotating Bose gases in the lowest 
Landau leveP^. Most recentljffi^, a generalised version of the composite fermion formalism, 
including a (pseudo)spin degree of freedoirP was applied to two-species Bose condensates, 
which can be realized experimentally in various wayd^^*^. Remarkably, in the boson studies, 
it turned out that the CF formalism produces close to exact wave functions even at the lowest 
angular momenta (typically smaller the the number of particles), i.e. far outside the actual 
quantum Hall regime for which it was originally intended. One particular issue that was 
noted for low angular momenta both in the single species casd^^®, for two-species bosond^, 
and in fact also in the context of highly excited states of electronic quantum Hall stated^^*^, 
is that the CF formalism frequently produces too many candidate states, or rather “hid¬ 
den” linear dependencies. More precisely, the number of seemingly independent CF Slater 
determinants before lowest Landau level (LLL) projection obeying the pertinent physical 
constraints [desired total angular momentum and (pseudo)spin quantum numbers, transla¬ 
tion invariance etc] is often considerably larger than the number of linearly independent CF 
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states after projection. For example, for 2+6 particles and angular momentum 4, there are 
41 distinct pairs of CF Slater determinants obeying all the physical constraints. However, 
after performing the projection - basically since various combinations of derivatives acting 
on the terms of the Jastrow factor, can conspire to result in the same polynomial - it turns 
out that there are only three linearly independent ones (which coincides with the actual 
dimension of the relevant eigenspace)^^. This is a mathematical feature of the CF formalism 
that remains to be fully understood, and the present paper aims to make progress in this 
direction. A systematic understanding of this issue is obviously of practical interest: numer¬ 
ical calculations would be signihcantly simplihed if one were able to systematically identify 
the set of linearly independent CF states a priori, from the unprojected form of the CF 
Slater determinants, rather than having to carry out the projection brute force, and looking 
for linear dependencies among the resulting polynomials. This is particularly important in 
the low angular momentum regime, where LLL projection amounts to a very large number 
of derivatives, and is thus computationally heavy. 

The main part of this paper concerns the systematic classihcation of linear dependencies 
in the CF description of two-species Bose gases with homogeneous interaction. We study 
the angular momentum regime L < N ■ M where N and M denote the particle numbers 
of the two species. In particular we focus on what we refer to as simple states. These 
were identihed in ReflTTI as a certain subset of all CF candidate states at a given number 
of particles and angular momentum. Technically, they are characterized by having at most 
one composite fermion occupying each A-level in the CF Slater determinants, and they can 
be shown to minimize the A-level “cyclotron” energy of the state. While diagonalisation 
within the full set of CF candidate states recovers basically the entire yrast spectrum for 
low angular momenta exactly or near exactly, the simple states still give a very accurate 
description of the low-lying part of the yrast spectrum, with typical overlaps > 97% for 12 
particled^. Even when restricting to simple states, the number of seemingly different CF 
candidates still tends to get vastly larger than the number of linearly independent states 
after projection. The reason for focussing on simple states hrst is that the concepts and 
techniques we use to systematically reveal linear dependencies are most easily introduced 
this way. However, many of these ideas also apply more generally to non-simple states 
(including single-species), as well as fermionic systems. These cases will be addressed later 
in the paper. 


3 


To summarize, our goal is to develop a systematic way of revealing linear dependencies 
at the level of the Slater determinants themselves, so that a minimal set of linearly in¬ 
dependent CF basis states is identihed before explicitly performing LLL projection. Our 
approach involves a letter string notation to represent the occupation patterns of the CF 
Slater determinants. Exploiting the fact that all so-called compact states (at some given 
angular momentum L) are translationally invariant leads to linear relations between states 
at L — 1. Further identities are found from a generalised version of translation invariance, 
as well as various ways of reordering the occupation patterns of the Slater determinants. 
All identities are proven analytically and supported by numerical calculations. Since this 
part of the paper is rather technical, we we will illustrate with examples, and defer some 
details of the proofs to an Appendix. We start by summarizing some necessary background 
theory in section [TT| Sec uni treats in detail the case of simple states, while general compact 
states and the fermionic case are discussed in sections W and |V| respectively. We end with 


summary and future perspectives in section VI 


II. TWO-COMPONENT ROTATING BOSE GASES 

We here give a quick summary of the model for two-species Bose gases in the lowest 
Landau level with homogeneous interaction, including their description in terms of composite 
fermions. For a more detailed introduction we refer the reader to ReflTTl The Hamiltonian 
for our system, two species of bosons in a two-dimensional harmonic trap of strength u, 
rotating at frequency H, is 

N+M / 1 \ N+M 

^ = Y. ( ^ + -^hj+ Y 27rc/(5(ri - r^). (1) 

i=l / i<j=l 

Here M denotes the number of particles of the majority species, and N is the number of 
particles of the minority species. The single-particle angular momenta are denoted by /*. 
We have assumed that all particles have the same mass m, and that the strength of the 
contact interaction, g, is independent of species. This is what we refer to as homogeneous 
interaction. In the weak interaction (dilute) limit this reduces to the well known lowest 
Landau level problem?^ in the effective magnetic held 2muj, 

N+M N+M 

H= {u - n)li + 2Tig 5{r]i-r]j). (2) 

2=1 2 <J =1 


4 




In the ideal limit (cu — f2) —)■ 0 one gets flat Landan levels, and so all the many-body physics 
of the system is determined by the interaction. Here rjj = xj + iyj are the dimensionless 
complex positions of the particles in nnits of the “magnetic” length hj{2mu). The second 
sum runs over both species since the interaction is homogeneous. Working in symmetric 
gauge, the lowest Landau single-body eigenstates with angular momentum I are 

'ipo^iiz) = Niz^ exp{-zz/4) l>0 (3) 

where we will suppress the Gaussian factor for simplicity from now on. A generic many-body 
wave function with fixed total angular momentum L is then a homogeneous polynomial of 
degree L, symmetric in the coordinates of each species separately. As previouslj^^^^ we will 
focus on translationally invariant states, i.e. polynomials invariant under a simultaneous, 
constant shift {K) of all coordinates, 

'^{z + K,w + K) ='^{z,w). (4) 

where z and w denote the sets of coordinates of the minority and majority species, respec¬ 
tively. 

Due to the species-independent interaction strength, the Hamiltonian is invariant under 
change of species, which implies a pseudospin-1/2 invariance. In Ref JTTl we explained how 
this spin analogy can be exploited to greatly simplify the analysis of the system’s many- 
body spectra. This aspect is of less importance in the present paper, where we will focus on 
the mathematical properties of CF states that are already known from Ref JTTl to be good 
candidates for the low-lying states of the system. 

A generic CF trial wave function for the bosonic two-species system is of the formpi 

^CF = 'PLLLi^Z^wJiz^w)) (5) 

where $ 2 , are Slater determinants for each species of the non-interacting CFs. They 
consist of the single-particle states 

= Nn,mZ"^L'^ , m > -U, (6) 

where is the associated Laguerre polynomial, and A„,m is a normalization factor. J is a 
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FIG. 1: Sketch of CF Landau (or “A”-) level occupancies for the pair of Slater 
determinants of a simple state. Such states are characterized by at most one composite 
fermion of a given species occupying any Lambda level. 

full Jastrow factor involving both species, 

V+M 

j{z,w) = n ivi-vj) 

i<j=l 

N M N,M 

= n n (yWk-wi) W{zi-wk). (7) 

i<j=l k<l=l i^k=l 

Vlll denotes projection to the lowest Landau level. The original, most direct projection 
method amounts to replacing the conjugate variables Zi, Wk by after moving them 

all the way to the left in the final polynomiaP. While other projection techniques are 
frequently employed in the literature in order to make numerical calculations less heavy, all 
our calculations are done using this original method. While the use of e.g. Jain-Kamilla 
projectioiP has been shown to make very little qualitative and quantitative difference in 
numerical calculations, our exact analytical results would presumably only be near-exact 
with this projection technique. We will briefly return to this point at the end of the paper. 


III. LINEAR DEPENDENCIES FOR SIMPLE STATES 

In this rather technical section we introduce, step by step, how to reveal linear depen¬ 
dencies between CF candidates at given particle numbers, total angular momentum and 
pseudospin quantum numbers. Each step introduces new linear relations between states, 
hnally leading to a systematic algorithm to reduce the set of CF candidates to a basis. This 
section focuses on “simple” states, i.e. states with at most one composite fermion of each 
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FIG. 2: The solid line shows the number of naively independent simple CF candidates as 
function of angular momentum for 4+8 particles. The dashed line shows the number of 
linearly independent wave functions in the lowest Landau level. 

species occupying a given A-level n, with minimal angular momentum (m = —n). The 
latter ensures translational invariance of these states. An example is given in FIG. The 
polynomial part of the corresponding single-particle eigenfunctions is r]n-n = (neglecting 
the normalization constant) which translates to rjn-n = after projection. 

The possibility of having, and necessity of understanding, dependencies between simple 
GF candidate wave functions is apparent from FIG. Indeed, for low enough angular mo¬ 
mentum (L < 15 in the case of 4+8 particles), the number of simple candidates exceeds the 
dimension of the LLL basis, meaning that they cannot possibly form a linearly independent 
set. Based on this we infer that there exist dependencies between the seemingly independent 
candidates at these low angular momenta, and indeed this is also the case for higher values 
of L, as can be checked by performing the projection to the LLL. Understanding these de¬ 
pendencies is therefore necessary in order to efficiently use the GF construction to study the 
two-component Bose system. 
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A. Letter string notation 


The simple CF wave functions form a very restricted set, which makes it possible to 
represent them compactly using strings of letters. For a state with N particles of type Z and 
M particles of type hF, the single particle states will be on the form (n, m) = [x, —x) with 
a;G{0, + M — 1}. Higher values of x are excluded as the power of the corresponding 

derivative, d^, would exceed the highest possible power of any variable in the Jastrow factor. 
Since no single-particle level can be occupied by more than one CF of a given species, there 
are four possibilities: A state (z, —i) can be occupied by a Z, a W, both a Z and a W, or 
neither a Z nor a W. These scenarios are denoted by “Z”, “IF”, “P” (“pair”) and “P” 
(“hole”), respectively. A many-body simple state is then represented by the corresponding 
string of N+M zero-indexed letters, ordered by increasing n. For example, the wavefunction 
for the state in FIG. (with the left figure representing Z-particles, and the right part 
corresponding to IF-particles) can be represented as 


T({zi},{u;i}) 



d?, 

d?, 

9?, 

Wl 

W2 

W3 

d}, 

d}, 

d}, 

Wl 

W2 


^Wi 

^W2 

CO 


J{z, w) = PWHW. 


( 8 ) 


The terms in the Jastrow factor contain exactly i variables of order at least N + M — i. A 
non-zero CF wave function must therefore have at most i differentiation operators of order 
N + M — i or higher. This implies that no suffix can contain more P’s than P’s in the letter 
string representation. Equivalently, no prefix can contain more P’s than P’s. 


B. Translation invariance 

Simple states are compact and thus translation invarianlpl, in the sense that their poly¬ 
nomial part satisfies 

T({^,}, {m,}) = T({^, + P}, {wi + K}). (9) 

We dehne the differentiation operators 

N N 

= ^dzi, Aiy = ^ dy,. . 
i=l j=l 
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( 10 ) 






Now, an equivalent way of stating translation invariance is 

(A^ + = 0. (11) 

These differentiation operators act on the CF polynomials by raising single particle states 
(i, —i) to (i + 1, —i — 1). For the polynomial PWHW, this gives 

(Az + Aw)PWHW = WPHW + PHWW = 0. (12) 

Applying translation invariance on states of angular momentum L thus gives linear depen¬ 
dence relations for states of angular momentum L — 1. Other examples are 

(Az + Aw)PPWZHH = PWPZHH + PPWHZH + PPHPHH = 0 (13) 

and 

(Az + Aw)PWZWH = WPZWH + PWHPH + PHPWH + PWZHW = 0. (14) 

These dependence relations reduce the number of candidate states, but not generally to 
a linearly independent basis set, except in some special cases when the total number of 
particles is low. 


C. Blocks 


A further reduction of the number of CF candidates, explored in this subsection, is due 
to invariance of the final CF polynomial under certain ways of permuting the occupation 
patterns of the Slater determinants. As mentioned, a string will represent the zero polyno¬ 
mial if it has a prefix with more FT’s than P’s and is therefore not interesting. An interesting 
situation occurs, however, when a prefix has an equal number of P’s and P’s. If the prefix 
is of minimal, positive length we call it a block. We repeat this with the remainder of the 
string until the whole string is partitioned into blocks. We write such a partitioned string 
as 


Xi o X 2 o ... o Xn, 


(15) 


where the Xi are blocks, and o denotes concatenation of strings. For example, the resulting 


polynomials in equation (12) have the blocks 


WPHW = WoPHoW, PHWW = PH oWoW. 


(16) 
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This partitioning of the string is useful because applying the Slater determinants to the 
Jastrow factor will only give non-zero terms if the differentiation operators from the hrst 
block are applied to the variables with the lowest exponents in the Jastrow factor. The 
differentiation operators from the second block are then applied to the lowest remaining 
exponents, and so on. This is used in Appendix|^to show that the corresponding polynomial 
is invariant under permutations of blocks up to a sign. The commutation rule is 


where the Xi are blocks, Ni denotes numbers of Z particles, and Mi numbers of W particles 


in block Xj. Incidentally, the linear dependence relation (|12|) can now alternatively be viewed 

( 18 ) 


as a consequence of the permutation rules for blocks, since 

WPH = W oPH = o kh = -PHW. 


Pure permutation invariance is not the only useful consequence of blocks; they also allow 
us to expand the concept of translation invariance. We note that, up to a combinatorial 
factor, the polynomial of a CF state with multiple blocks equals the symmetrization of the 
polynomials of the blocks when they are considered as individual states. Since each block 
is translation invariant, we can differentiate the blocks separately and thus obtain more 
dependence relations. For example, since 


{Az + Aw)PWZH = {Az + Aw)PWH = 0, (19) 

the following dependence relations can be derived from the state PWZPfPWPf : 

{PWHZ + PH PH) o PWH = PWHZPWH + PHPHPWH = 0, (20) 

PWZH o {WPH + PHW) = PWZHWPH + PWZHPHW = 0. (21) 

In contrast, using translation invariance of the whole state would merely give 

PWHZPWH + PHPHPWH + PWZHWPH + PWZHPHW = 0. (22) 

In general, translation invariance for individual blocks means that if we have a dependence 
relation 


= ( 23 ) 

i 

where the Xj are letter string representations of wave functions, then it is also true that 

E^aOWoXfc = 0, (24) 

i 

where X^, X^, are arbitrary letter string representations of wave functions. 
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D. Reflection 


In this section we prove that CF polynomials of simple states are invariant nnder what we 
refer to as reflection symmetry. This symmetry is somewhat related to translation invariance, 
bnt it is not captured in the dependence relations found above. Reflection is defined as 
follows. Assume X is a letter string representing a CF wave function or block. Then we 
define the reflected string X' as the result of the following operations: 

1. Invert X, i.e. ABCD —)• DCBA. 

2. Swap H -H- P. 

Some examples are 

PWPZHH' = PPZHWH, PPWZPHZWHH' = PPWZPHZWHH. (25) 

Counting from the beginning of the reflected state, the Z(W) particles occupy the positions 
that were vacant of W(Z) particles in the original state when counting backwards from the 
end. It is not difficult to see that reflection leaves X, M and L invariant, and that X" = X. 
The non-trivial result is that 

X' = (-1)^^+^X. (26) 

We need a small lemma for the proof: 


{AzXy = AwiXy (27) 

Applying Az and then reflection to X can be described by the list of operations 

1. Create the sum over every possible Z moving to the right. 

2. Invert X, i.e. ABCD —)• DCBA. 

3. Swap H -H- P. 

Or, equivalently 
1. Invert X. 
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2. Create the sum over every possible Z moving to the left. 

3. Swap H -H- P, or equivalently, (missing Z) —>■ W, (missing W) —)■ Z. 

Further, since a Z moving to the left is the same as a missing Z moving to the right, we get 
that this procedure is also equivalent to 


1. Invert X. 

2. Swap H ^ P. 

3. Create the sum over every possible W moving to the right. 


and this hnal list is just the procedure that describes Aw{X'). 


We are now ready to prove reflection symmetry (26) by induction on L\ 


For L = 0, every possible string representation will consist of N Z’s and M IF’s in some 
order. Since there are no P’s or P’s, the reflected string will just be the inverse string. We 
know from section IIII Cl that W and Z blocks anticommute while Z commutes with Z and 
W commutes with W. Since one possible way of commuting the blocks so that X —)■ X' 
involves moving every block past every other block, we get a sign of exactly (—1)^^ from 
interchanging neighbouring hF’s and Z’s. This shows that the result is satished for L = 0. 


Now, assume that (26) holds for L — 1 and consider the following expression. 


Az{X - (-1)^^+^X') = AzX - (-l)^^+^Az(X') 

= AzX + {-l)^^+^Aw{X') = (AzX) - {-l)^^+^^-^\AzXy = 0. (28) 


The second equality is due to translation invariance, A^X = —AyyX, then we use lemma 


(27) and hnally the induction hypothesis for L — 1, since A^X is a sum of strings with 


angular momenta L — 1. Since Az = —Aw, we also have 


Aw{X — (— 1 ) 


NM-\-L \^/ 


X') = 0. 


(29) 


Since X-(-l)^^+^X' is a symmetric function in both z variables and w variables, Eqs. 
(28 29) imply that X — (—l)^^+^X' is constant in all its variables. Since L > 0 it must 


vanish. Thus, (26) holds for L, which completes the proof. 
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E. Generalised translation invariance 


Applying the results of block permutation, reflection symmetry and translation invariance 
does reduce the number of CF candidate states signihcantly, but still leaves dependencies 
that cannot be explained from these symmetries. The hnal piece to the puzzle is a generalised 
version of translation invariance. We dehne differentiation operators 


N Oa M Oct 

= Awa = J2 




a e N. 


(30) 


i=i j=i 

These operators commute with the Slater determinants, since the latter consist only 
of derivatives for simple states. They can therefore be applied directly onto the Jastrow 
factor in order to study their action on the simple states. Since the Jastrow factor is anti¬ 
symmetric in all its variables, and Azn+Awn is symmetric in the same variables, the product 
{Azn + A]Vn)J{z,w) is auti-symmetric too. However, since the Jastrow factor contains all 
powers from 0 to iV J- M — 1, any non-zero resulting term from the differentiation must have 
two exponents that are equal. The anti-symmetry then ensures that the result is zero. 

The simple states that result from this higher order differentiation can be obtained in 
the same way as for translation invariance, using letter string notation. Applying Aza to a 
string results in a sum over every possible move of a Z a places to the right. A difference 
from translation invariance is that such a move can cause a Z to move past another Z, and 
this causes a minus sign. The reason is that the polynomials that the strings represent are 
dehned with differentiation powers occurring in increasing order in the Slater determinants. 
For example, we have 


(Az2 + Aw2)PPWWZHH = -WPPWZHH + PWWPZHH - PPHWPHH + PPWHZWH 

= 0. (31) 

As with regular translation invariance, (Aza + Awa) may also be applied to individual 
blocks to create one dependence equation per block in a state. In addition to producing 
more dependence equations, this approach has the added beneht of constraining the values 
of a that need consideration. To see this, consider a string composed of blocks 


Xi o . . . O Xj o Xj +1 o ...oXn. (32) 

If the length of the block Xi is /*, the application of, say, Aza on this block will necessarily 
move a Z from X^ to Xj+i if a > li. Such a move will either turn a P into a W, or a Z into 
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an H in Xj. Bnt then, the prehx Xi o... oXj contains more if’s than P’s, and conseqnently, 
the string represents the zero polynomial. The conclnsion is that any term where a Z or W 
moves from one block to another in the original string resnlts in a vanishing polynomial. In 
particnlar, when acting with Aza + Aw a on a block Xi with length /j, only a < U will give 
something non-zero. 


F. Combined results 


We are now in a position to describe the main resnlts of this work. In this section we 
present an algorithm that rednces the nnmber of CF candidates for a given N, M, L by 
as mnch as possible, according to the known linear dependence relations detailed above. 
The algorithm is easy to implement becanse the set of all simple states is closed nnder the 
operations leading to dependence relations. 

The algorithm can be snmmarized in the following steps: 


1. Constrnct all simple CF candidates at N,M,L. This set of candidates is named 
{Tcf}- 

2. Constrnct dependence relations by applying all combinations of permntation and re¬ 
flection of blocks. Use the dependence relations to rednce the set {Tci?}. 

3. For all AL = 1,..., N + M — 1, do: 


(a) 

(b) 


Constrnct all simple CF candidates at N, M, L + AL. 

Constrnct additional dependence relations nsing the generalized translation in¬ 


variance condition (39) with a = AL, applied to all the blocks of length greater 
than a of the states in step 3(a). Simplify the new dependence relations using 
those found in step 2. 


4. Use the relations found in step 3 to further reduce the set 


Using the algorithm presented above, we can compare the size of the reduced {Tcf} to 
the actual number of linearly independent wave functions, found by a brute-force calcula¬ 
tion. We hnd that {Tcf} is fully reduced in all cases, i.e. that the reduced {Tcf} is a 
linearly independent set of wave functions, and thus constitutes a basis for the simple CF 
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candidates. The algorithm thus succeeds in removing all linear dependencies without doing 
any projections to the LLL, i.e. without actually comparing the very complicated polyno¬ 
mials in the final wave functions. The statements above have been verified numerically for 
up to 14 particles. 


To illustrate the benefit of applying the algorithm above, FIG. 3a shows the number of 


naively independent, and actually independent, simple CF candidates, while FIG. [3b] show 
their ratio, for 4-|-8 particles. For low L the size of the basis is two to three orders of 
magnitude smaller than the size of the naive set of candidates, meaning that projecting all 
of them would be extremely wasteful, computationally. For somewhat larger L, the ratio 
decays rapidly with L, but the absolute number of states is still significantly decreased by 
the algorithm. 

We can see how reflection symmetry, block permutation, ordinary and generalised trans¬ 
lation invariance contribute to the explanation of linear dependencies. FIG. shows the 
number of candidates after removing dependencies caused by the different mechanisms. 
First, we notice that overall, reflection seems to contribute the least, particularly for small 
L. This is easy to explain: reflection removes one out of two blocks that are reflections 
of each other, but only for blocks that are not reflected unto themselves [PPWHH is an 
example of a block equal to its reflection). For low L, many blocks are short (contributing 
no or little angular momentum to the total), and given that the shortest non-self-refiected 
blocks are the four-particle blocks PZWH and PWZH, reflection provides few dependence 
relations. Second, block permutation is increasingly important with decreasing L for the 
same reason: as L decreases we can form more and more small blocks, increasing the number 
of candidates that only differ by the order of the blocks. On the other hand, the largest L 
for which one can have more than one block is L = N{M — 1). This is the two-block state 


P...P W...W H...HoW 


(33) 


N M-N-l N 

For L > N{M — 1), all candidates are one-block states, and permutation symmetry tells us 
nothing. 

Last, we notice that the number of candidates after generalised translation invariance is 
symmetric about Lmid = NM/2. In fact, the generalised Tl-line equals the fully reduced 
line from FIG. and this is also the case for the other numbers of particles we have 
studied numerically. However, remembering that reflection and block permutation can be 
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N=4,M=8 



N=4,M-8 



FIG. 3: (a) Logarithmic plot of the number of naively independent, and actually 
independent, simple CF candidates, for 4+8 particles, (b) The corresponding ratio of the 
number of naively independent to linearly independent simple CF candidates. 
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N=4,M=8 



FIG. 4: Logarithmic plot showing the number of seemingly independent simple states after 
reducing the naive set by the mechanisms described in this section, again for 4+8 particles. 

TI stands for translation invariance. 


exploited very easily before applying the generalised translation invariance conditions, this 
is a valuable, and often necessary, thing to do: especially for small L, neglecting to exploit 
block permutation will lead to extremely many unnecessary generalised Tl-equations. The 
symmetry about Lmid lends hope to the idea that there might exist an analytic combinatoric 
formula for the dimension of the simple CF basis at any N, M, L. We comment on this in 
section m 


IV. GENERAL COMPACT CF STATES 

The simple state with the highest possible angular momentum L for given N, M is pro¬ 
duced by hlling the “ladder” of CF orbitals {(n,—n)} from below, giving the letter string 
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representation: 


P_^ W. . .W, (34) 

N M-N N 

This state has angular momentum L = N ■ M. For L > N ■ M, no simple states exist. The 
relevant CF candidates for the low-lying eigenstates with angular momentum L > N ■ M 
are known as compact stated. A single species CF candidate is said to be compact if its 
Slater determinant is such that: 

1. for any occupied CF orbital t/’n.m, m > —n, the orbital is also occupied, and 

2. for any A-level n occupied by k CFs, the level n -|- 1 is at most occupied hy k + 1 CFs. 

As mentioned earlier on, simple states satisfy these requirements trivially, i.e. simple 
states are special cases of compact states. Compact wave functions have some well-known 
propertied^. First, compact states obey translation invariance of the polynomial part of 
the wave function, like simple states do. Second, one can use row reduction on the Slater 
determinants after projection to show that the effective CF orbitals can be writteiP 




n,m 


.Zi = z, 


m-\-nQn 


(35) 


For a two-component compact CF candidate, the Slater determinants of each species are 
required to be compact separately. The compact state with the largest obtainable L has all 
particles sitting compactly in the lowest A-level in both Slater determinants, so this state 
is nothing but the “221” state, ^/> 22 i = Ili<j{zi - Zjfl\k<i{wk - wiyilm<n{zm - Wn). The 
angular momentum of this state is L 221 = N{N — 1) -|- M(M — 1) -|- A^ ■ M. The compact 
states are thus our objects of interest in the angular momentum range N ■ M < L < L 22 i- 
We now consider the generalizations of the linear dependence concepts presented for 
simple states to general compact states. Apart from certain modihcations which will be 
described, all the concepts developed for simple states will also apply to compact states, 
except for reflection symmetry, for which we have not found a suitable generalization. The 
modihcations are necessary due to the fact that the set of compact states is not closed 
under the action of (generalised) translation invariance or permutation of blocks. These 
modihcations will be the topic of this section. 
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A. Notation 


The letter string notation is not adequate to describe compact states. In this section, we 
will, in addition to simply listing which CF orbitals are occupied or writing down the Slater 
determinants, employ a diagrammatic representation as follows. The letters P, Z and W will 
signify occupation by both species, only the Z species, and only the W species respectively, 
and the possible orbitals are shown in a m vs. n diagram. Unoccupied orbitals will be 
displayed as dashes. For instance, the CF candidate which has '^ 2,-2 and '^ 3-3 occupied by 
the Z species, and 'ipofl, 'ipo,!, and iIja-a occupied by the W species will be represented 
by the A-level diagram 


5------- 

4 w - - - - - 

3 z - - - - 

n 2 P — — — 

0 W W 

-5 -4 -3 -2-10 1 

m 


(36) 


All allowed A-levels n = 0 ,...,A + M — 1 are included in the diagrams for completeness, 
making it easier to compare candidate states, but they may be truncated after the largest 
value of m. 


B. Generalised translation invariance 


We wish to exploit that the Jastrow factor vanishes under the generalised translation 
operator: 


i^Za + ^Wa) J 


N 


M 


Y.K + Y.fK., J 


\ i=l 


i=i 


0 V a e N 


(37) 
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The product of one projected Slater determinant, say for the Z species, and the Z-part of 
the generalised translation operator, is 


k=l 

N 


'>, 1=1 


= E(-ir'n 




N 


"'2fe 


T.31. 


k=l 


Vi=l 


N 


N 


k^i 


(38) 


i=l \ cr 


N 


N 


~E(E( V'n^(i)+a,m^(i)-a(^i) n t/’' 


^(T(fc) >^Cj(fc) 


Ut) 


i=l \ cr kf^i 

where a is a permutation of N elements, and (—is the signature of the permutation. 
Thus we see that ) is equivalent to the sum of replacements of a single 

orbital 'ijjn,m —t 'ipn+a,m-a- In the diagrammatic notation, this amounts to the sum of all 
moves taking a Z (a W) a slots “up and to the left” in the A-level diagram. If we take to 
mean the Slater determinant <h where the orbital '0n,m in row i of the determinant matrix 
has been replaced by 'ilJn+a,m-a, then 

N M \ 

+ '^i^z){^w)‘j J = ^Z^W i^Za + ^Wa) J = 0 (39) 

>i=l j=l / 


due to (37). This means that we can use generalised translation invariance to find dependence 


relations as we did for simple states. Notice that (39) is not equivalent to having Aza + ^Wa 


to the left of the Slater determinants. For compact states we generally do not have 

{Aza + Awa) ^Z^wJ = 0 O > 1. 


(40) 


re- 


We will now see that some care needs to be taken when applying this, since some 
suiting terms may not be compact. Consider the one-component CF candidate with Slater 
determinant 

(w?a° J (w^^o, 

— 

i.e. two CFs in the lowest A-level. The action of this determinant on Aw is 


(41) 


^wAw = 


J {wldl 


+ 




(42) 
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and since the two resnlting determinants are not compact, they are not even CF Slater 
determinants: the CF orbitals for a non-compact state have the general form ([^, which these 
determinants clearly do not. However, this will not prevent ns from applying translation 
invariance. We will simply inclnde these non-CF contribntions as nnknown fnnctions Xi 
the same footing as the compact states when writing down the linear relations, and attempt 
to eliminate them when performing the row rednction to independent wave functions. In the 
worst case scenario, all equations will still contain these non-CF states after row reduction, 
and we will have learned nothing about dependencies within the set of compact candidates. 


However, the results presented in IVD indicate that, at least for the lowest sub-band of 


compact CF states, this never occurs, i.e. we are in fact able to hnd all dependencies, and 
thus reduce the set of candidates to a proper basis, as was the case for simple states. 

C. Blocks 


From the form of Eq. (35), it follows that the appearance of non-zero powers of the 


coordinates Zi, wj in the CF orbitals does not affect which terms in the Jastrow factor that 
survive the differentiation imposed by the Slater determinants. Therefore, the classihcation 
of blocks is also possible for compact states, and is simply achieved by counting derivatives 
as for simple states. For instance, the Slater determinant pair represented by the diagram 


in (36) splits into two blocks: the hrst block Xi contains A-levels 0 and 1, and the second 


X 2 contains A-levels 2 through 5. Notice that X 2 is in fact a simple block, and that we may 
apply reflection symmetry to X 2 to immediately acquire the relation 

/ . \ / , \ 


5 
4 
3 

n 2 
1 
0 


W 


Z - - - - 

P - - - 


\ 


-5 -4 -3 -2 
m 


W W 
-1 0 1 


= (-!)■ 


/ 


5 
4 
3 

n 2 
1 
0 


W - - - - 

P - - - 


V 


-5 -4 -3 -2 

m 


W W 
-1 0 1 


(43) 


/ 


In general, we may apply reflection symmetry to all simple blocks of a candidate, even 
though the full state is not simple. This is as far as we have been able to pursue the 
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reflection concept for compact states. 

Now we generalise permutations of blocks to compact states. It should be clear that 
any block spanning more than one A-level (more than a single Z or W) necessarily ends 
on an empty level. However, a block may begin on a level containing multiple occupancies, 
and in particular, the lowest A-level can contain any number of particles. Permuting two 
blocks may then violate the second condition for compactness. If this is the case, then the 
permutation of this block with any other block will result in a non-CF wave function, with 
elements not corresponding to the correct expressions for the orbitals. However, we may use 
this fact to get rid of some of the variables Xi that result from the translation invariance 
equations. To illustrate this, consider the action of the block 


4 _ _ _ _ _ _ 

3 Z - - - 

n 2 W — — — 

0 F W 

-4 -3 -2-10 1 

m 

on Aw 2 - Among other terms, this will produce 


(44) 


4 _ _ _ -- 

3 Z - - - 

n 2 W W - 

1 - (45) 

0 F - 

-4 -3 -2-10 1 
m 

but with the second W in A-level 2 given by Wid^., which is not the correct expression for 
the '^ 2,-1 orbital in this non-compact case. However, we realize that permuting the two 
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blocks of the compact candidate 


4 _____ _ 

3 p _ _ _ _ 

n 2 — — — — 

1 Z - - (46) 

0 W W 

-4 -3 -2-10 1 


m 


will give exactly (45) apart from a minus sign given by Eq. (0. We conclude that some 
terms resulting from applying Azn or Awn may at hrst not seem compact, but turn out to 
be just that by carefully considering permutations of blocks. It is this fact that allows us to 
eliminate many unknown functions Xi from the set of dependence relations, which in turn 
gives us more information about the dependencies between the compact candidates. 


D. Combined results 


The modihcations described above are reflected in modihcations to the algorithm de¬ 


scribed in section III F To obtain results for the lowest part of the yrast spectrum, we will 
consider the CF candidates that, at a given N,M,L, minimize the CF cyclotron energy 


N+M 


(47) 


2=1 

where rij are the A-level quantum numbers of the occupied CF orbitals. We immediately 
notice that the application of Slater determinants on Aza + ^wa raises K and lowers L by 
a units, and that permutation of blocks and reflection of simple blocks leave K invariant. 
Remember also that the action of a block with N + M particles on -|- A^y^ trivially 
vanishes if a > A -|- M, because too many derivatives will act on the Jastrow factor. 

The algorithm can be summarized in the following steps: 


1. Construct all compact CF candidates at A, M, L with minimal CF cyclotron energy, 
Kmin- This set of candidates is named 

2. Construct dependence relations by applying all combinations of permutation and re¬ 
flection of blocks (reflection only for simple blocks). Use the dependence relations to 
reduce the set 
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3. For all AL = 1,..., iV + M — 1, do: 


(a) Construct all compact CF candidates at N,M,L + AL with CF cyclotron energy 


^ ^min AL. 


(b) Construct additional dependence relations using the generalized translation in¬ 


variance condition (39) with a = AL, applied to all the blocks of the states 


in step 3(a). Simplify the new dependence relations using those found in step 
2. Generally, some dependence relations will simplify to 0 = 0, and some will 
contain unknown functions Xi- 


4. Using Gaussian elimination or similar methods, eliminate as many of the Xi as possible. 

5. Use the dependence relations not containing any Xi fo reduce even further. 

Again we compare with a brute-force calculation of the number of linearly independent 
states. Surprisingly, we hnd that is fully reduced in all cases, i.e. that the reduced 

is a linearly independent set of wave functions, and thus constitutes a basis for 
the GF candidates in the lowest sub-band K = Kmin- We have checked this for all 2072 
combinations of N, M, L where iV -|- M = 2,..., 12; iV = 0,..., [{N -|- M)/2J ; L = N ■ M + 
1,..., F 221 . An example of the numbers of candidates before and after exploiting dependence 
relations is seen in FIG. for 3-1-7 particles. A clear pattern is visible, where the two lines 
agree in the cusps at certain values of L = L^. The values of Lj are characterized by having 
a unique state with minimal GF cyclotron energy Ki such that, at L = Lj -|- 1, there exist 
at least one state with the same K = Ki. This happens when it is not possible to reduce K 
by increasing L\ at L not equal to any of the L*, K strictly decreases with increasing L. 

For higher sub-bands K > Kmin, we hnd some, but not all, linear dependencies. Specih- 
cally, there exist dependencies between linear combinations of states in one band, and linear 
combinations of states in lower bands. As our methods only produce dependence relations 
between states with equal K, we do not capture inter-band relations. We do, however, 
capture all intra-band dependencies. We expect that one would need to explicitly treat the 
non-compact contributions Xi order to fully understand inter-band dependencies. This is 
left for future work on the subject. 
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N=3,M=7 



FIG. 5: The number of compact candidates before and after exploiting dependence 
relations, for N = 3, M = 7. The dashed line shows the naively independent states, while 
the solid line shows the remaining candidates after removing dependencies. For 
42 < L < L 221 = 69 the two lines agree, i.e. the candidates are linearly independent to 

begin with. 

V. ELECTRONIC CF STATES 

We briefly comment on the application of this work’s hndings to CF states for the standard 
quantum Hall effect - one-component electron systems in magnetic fields. We maintain the 
disk geometry and the projection method used so far, and comment on other geometries and 
projection methods later. The CF candidate wave functions in the disk are given by 

= Vlll ($ (48) 

with the usual interpretation that each electron captures an even number, 2p, of vortices to 
form composite fermions. In contrast to the bosonic case, an even power of Jastrow factors, 
.P in the simplest case, does not only contain terms with one set of exponents like does. 
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Therefore, there does not seem to be a straightforward way to identify blocks in a manner 
similar to what we have presented in this work. The largest allowed derivative in the Slater 
determinant is also higher now, so even if only the m = —n CF orbitals are occupied, there 
will be more H^s than P’s in the string letter description, and we have not been able to give 
meaning to the reflection operation in this case. 

For bosons, the result of acting on a (compact) CF candidate with the translation operator 
Az could be written as a sum of modihed determinants acting on the unmodihed Jastrow 
factor because J itself is translation invariant. This is no less true for and since the 
form of the Slater determinants is the same as before, translation invariance will reveal linear 
dependencies also for electronic CF wave functions, as long as they are compact. However, 
generalised translation invariance does not hold for 2p > 1: 

AzaJ^^ ^0 a > 1, 2p > 1 (49) 

To see how Az acts on a Slater determinant, we write 

AzVlll ($ = Vlll ( 50 ) 

Before projection, 2 : and ^ commute, so we can evaluate the unprojected part of the right 
hand side of the equation. We need the following property of associated Laguerre polyno¬ 
mials: 


= (n + m)L” (x) - (n + (x) 


( 51 ) 


From this we easily hnd 


(52) 


ZllJn,m{z,z) = ZL'^{ZZ)Z'^ 

= z^-hzL^{zz) 

= {n + m)L^-\zz)z^-^ - {n + l)L^-^\zz)z^-^ 

= {n + m)iJn,m-l{z, z) - {n+ l)'lpn+l,m-l{z, z) 

where we have ignored normalization of 'ipn^m and a factor 1/2 from the argument of L™. 
The interpretation is that AzVlll is equal to a linear combination of CF candidates, 

where each term has one CF moved either to the left, or up and to the left in the A-level 
diagram. Now, if <h is compact, then 


AzVlll ($ J '"’) = 0 


(53) 
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We may use this to find linear dependencies between compact and non-compact states. For 
instance, the v = 1/3 state for iV = 3 has the CF-orbitals {n,m) = (0, 0), (0,1), (0, 2) 
occupied, and applying Az to this state gives 

Vlll ((-$1 - $2 - <h3) j^) = 0 (54) 

where $i,<F 2 ,‘h 3 are the three possible composite fermion particle-hole excitations with 
AK = 1, AL = -1. 

VI. DISCUSSION AND OUTLOOK 

The main result of this paper is a comprehensive algorithm that reduces the set of seem¬ 
ingly different, low-lying compact composite fermion candidates at given N, M, L to a basis 
set of linearly independent states, prior to LLL projection. This may amount to reducing the 
number of states by as much as orders of magnitude. The algorithm is based on exploiting 
translation invariance, and invariance of the final wave function under certain permutations 
of the Slater determinant occupation patterns, to identify the hidden linear dependencies. 
We presented the derivation in the context of so-called simple states for two-species bosons, 
and then outlined the pertinent modifications for general compact, one- and two-species 
bosons. In all cases, we found complete reduction to a basis when the lowest sub-band was 
considered. 

At this time we have not yet been able to fully prove that generalised translation invari¬ 
ance (possibly helped by reflection and block permutation) is sufficient to produce a basis 
of simple candidates for general N, M, L. However we have found a recursive formula for 
the number of linearly independent stated, only involving N, M, L, and work is underway 
to complete the proof that generalised translation invariance explains all the dependencies 
between the simple CF candidates. We hope that this will lead to a revised algorithm that 
directly produces bases for the space spanned by the simple states. 

We also commented on the case of fermions which is more complicated, since many of our 
results mathematically rely on the single-flux attachment of the boson case. A possible way 
around this might be to study fermionic CF-type wave functions of the form ^|JF = J ■ 'ipsi 
where J is a Jastrow factor. In other words, instead of the usual double flux attachment 
for fermions, leave one of the Jastrow factors “outside” the projection. Essentially this 
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would amount to a somewhat different LLL projection method, and one would have to test 
numerically if it produces good wave functions. This is left for future study. 

Although different choices of projection tend not to matter much when it comes to over¬ 
laps, energies etc. of CF trial states, they do produce somewhat different wave functions. 
There is thus reason to believe that using, e.g. Jain-Kamilla projection instead of “brute 
force” projection, the analytically exact linear dependencies of this paper might only be 
near-exact. However, it is possible that exact identities similar to those presented in this 
paper exist also for other projection methods, and that the ideas presented here may aid 
in identifying them. This, too, is left for future study. Similarly, it would be of interest to 
re-examine the issue in different geometries, in particular on the sphere. 

Another obvious question is whether there is some more direct, qualitative or intuitive 
way of predicting the linear dependencies from the CF Slater determinants - complementary 
to our rather mathematical algorithmic approach. We have not come up with a good answer 
so far. Indeed, some of the dependencies are so mathematically subtle that they may not 
have a simple, qualitative explanation within the CF phenomenology. 
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Appendix A: Proof of block permutation invariance 

This appendix details the proof that permuting two blocks in a CF polynomial will leave 
the polynomial invariant up to a sign as stated earlier (HIl- We consider a string with K 
blocks. The j’th block has uj Z’s and pj hF’s. The total number of Z{W) is N{M) as 
always. Define the variables 


Xi = ZiW 0 < i < N, Xi = Wi-N yN<i<N + M 


(Al) 


and exponents 


{tti I 0<i<iV},{ai 
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N <i< N + M} 


(A2) 


such that the first (second) set contains the exponents of the z{w) differentiation operators 
in increasing order. We can now use the symmetric group of n elements, Sn to write the CF 
polynomial as 

N+M N+M 

4- = E E (- 1 )“''*' n ^ n 4:'- m 

o'GSn^Sm '^—1 k=l 

The next step is to factorize each permutation into one permutation distributing the vari¬ 
ables in the different blocks and one permutation permuting within the blocks. We dehne 

s. = 0 , S,. = 0 S,.,, = 0 (A4) 

j = l j=l j=l 

and the quotient (not a group) 

S/ = Sn/S^ © Sm/S^, = { min({ a o p | a 0 S',, © S',, }) | p 0 S'at © S'm } , (AS) 

where min() takes the lexicographically smallest element. In words, for each distributition of 
{ 1, ..., iV + M } into 2K parts of lengths (z/i,..., Pi, ■■■, S/ contains one permutation 
with this distribution and the permutation is sorted in ascending order within each part. 
We can now factorize the determinant permutations as 

N+M fjai N+M 

= E E E (- 1 )"'+'''"''' n ^ n (A6) 

cr&S^ p&Siy+S/j. *=1 k=l 

We use that r 0 Sn+m must distribute Uj 2 :’s and pj re’s in the j’th block for the contribution 
to be non-zero. Let us dehne the permutation 

a' = (1,..., z/i, iV + 1, ..., iV + pi, z/i + 1,. . . , M) (A7) 

which is the lexicographically smallest such permutation. One then arrives at the factor¬ 
ization 

T = p o a' o a, p0 S'„+^, a 0 S/. (A8) 

Also, since we require that the indices distributed in the differentiation blocks are the same 
as the indices distributed in the variables, the a permutations must be the same. This allows 
us to write 

N-\-M N-\-M 

s = E E E n ^ n (A9) 

u&S/ p&Sv+S^ i=l '^(pocr)i k=l 
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Some further notation is needed to reduce this to an explicitly block order independent 
expression. We introduce the subgroups 

.7-1 K 

{Sa)j = ^ © Saj © ^ Ifii C Sa, Cl G { Z/, /i, Z/ + /i } , j G { 1, ..., } . (AlO) 

2 = 1 2=J'-|-1 

Also, for j G { 1,..., } we dehne 


Tlji 


V 0 < i < z/j 


i + ELl 

N + i - Vj + Eh=i y i^j < i < Vj + /Cj 


(All) 


such that riji is the i’th index of the j’th block. We can now write 


K 




^ 1—r —1 




X 


{Toa'oa)n.f, 


(A12) 


s = (-I)'"'' y n E i; (-i)|<'|+m n „., 

creS/i=ipe(s^)je(S^)j-Te(s„+^)j i=i k=i 

All the differentiation operators in the j’th block are of order at least Uji. We can carry 
out this differentiation in each block to obtain 


K ’^3+k-i '^7+1*7 

^ = (-l)KI ^ J]^(_l)|p|+M II ^ tt .fe-i 


n •^(TO(T'oa)n (U _ l)] 

aeS/3 = lP,r i=l k=l ^ 

The product over j can be carried out for the numerical factor and we get 

K 


X, 


{Ujk - 1)! 


. (A13) 


t = (-i)|-'|ce n E E (-1) 

a&S/ J = 1 p&{Sv)j®{Sii)j T£{Si,+fj,)j 


‘'j+H PtO-n-.-riji ^j+Pj 

|p| + |r| TT TT ^k -1 


n 


n 


...... „ a„..-nji J.J. ^(rocr'oo-), 

- ‘-i 


(M4) 


'jk 


where 


C 


nr=o 


N+M-l -I 


Z! 


U^iKL 


Uj+Pi-l 


%\ 


(A15) 


This factor is independent of the order of blocks, and the same is true for the reduced 
differentiation exponents, an^i — nji. The only block order dependent part of the polynomial 
is the sign (—1)I'^'L Finally, let us deduce the sign rule. 

It is sufficient to consider a wave function with only two blocks Xi 0 X 2 . The two different 
orders of blocks will have: 


a'l : Zi...znWi...wm -t Zi...z^^Wi...w^^z^^j,i...znw^^+i...wm 


(A16) 


and 


a '2 : Zi...ZnWi...Wm -t Zi...Zy^Wi...W^^Zy^+i...ZNW^^+i...WM 


(A17) 
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respectively. The signatures are 




The relative signature of the two a' is the product of the signatures: 


- 1)1 




which was stated in Eq. lUl- 


(A18) 


(A19) 
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